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Abstract 

For a simple model of mutually interacting qubits it is shown how the er- 
rors induced by mutual interactions can be eliminated using concatenated 
coding. The model is solved exactly for arbitrary interaction strength, for 
two well-known codes, one and two levels deep: this allows one to see under 
which circumstances error amplitudes add coherently or incoherently. For 
deeper concatenation, approximate results are derived which make it possible 
to calculate an approximate "threshold" value for the product of interaction 
strength and free evolution time, below which the failure probability for an 
encoded qubit decreases exponentially with the depth of the encoding. The 
results suggest that concatenated coding could fully handle the errors arising 
from mutual interactions, at no extra cost, in terms of resources needed, from 
what would be required to deal with random environmental errors. 
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I. INTRODUCTION AND MODEL 



In previous papers |TJJ2j I have discussed how unwanted mutual interactions between 
qubits can be a source of error in a quantum computer (or quantum memory register) and 
studied some of the characteristic features of this type of error. It was found in |l| that, if 
the characteristic interaction energy is hS, the fidelity T of the register decreases with time 
t and with the total number of qubits N as 

T ~ e~^ 2 (1) 

Equation (|]) suggests that (5t) 2 can be interpreted as a sort of failure probability per qubit. 
In it was shown that the general scaling ([!]) with t holds even when the register is being 
acted on by "gates" such as one might use in the course of a quantum computation: that 
is, the error accumulates quasi-coherently in time, and there is no "quantum Zeno effect" 
associated with normal gate operation. (This suggests, in particular, that it is not necessary 
to simulate a full computation to estimate the magnitude of this effect; instead, it should 
be enough to calculate the survival probability of an appropriate initial state of the register 
under the action of the appropriate Hamiltonian [[J.) 

Mutual interactions are Hamiltonian and hence, in principle, reversible, and there exists, 
in fact, a variety of special techniques for eliminating, or greatly reducing, these errors: in 
particular, techniques similar to, and inspired by, the use of refocusing pulses in magnetic 
resonance 0|| . Nonetheless, it is worthwhile to look at the question of how well these errors 
could be eliminated by the use of ordinary error-correction methods based on quantum error- 
correcting codes H[7|]! and in particular on concatenated coding, if only because one would 
like to be able to make an informed choice between various available methods. Moreover, 
although pulse methods could in principle be used to deal with certain kinds of environmental 
noise, it seems very likely that, in large-scale implementations of quantum computers, error- 
correcting codes will have to be used, in any event, to deal with at least some of these 
random environmental errors. Assuming, then, that these codes would already be in place, 
it is natural to ask how much one could get from them. 

While there is no question that the present theory of quantum error-correction codes 
covers mutually interacting qubits as a special case, this is not a case that is often explicitly 
discussed, and, in particular, it is not immediately clear whether, in order to obtain a 
threshold estimate for these kinds of errors one should add them coherently (adding error 
amplitudes) or incoherently (adding error probabilities). This is one of the issues explored in 
this paper, for the two best-known quantum error-correcting codes: the 5-qubit and 7-qubit 
codes, and their concatenated forms. The conclusion, at least for the kind of interaction 
Hamiltonian assumed, is that these codes could fully handle the errors arising from mutual 
interactions, almost certainly at no extra cost in terms of resources needed (e.g., depth of 
concatenation). This is an encouraging result which may not have been immediately obvious 
before. 

The interaction to be considered here is one that results in pure phase shifts, conditioned 
on the relative state of neighboring qubits, according to a Hamiltonian 

N-l 

H =Y, h$ a zn°z,n+l (2) 
n=l 
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Examples of physical systems leading to interactions of this form have been discussed in 
previous publications in @] it has also been shown that, for instance, changing a z to 
a x (bit errors instead or phase errors) does not modify substantially the scaling ([1]). 

In order to simplify the calculation, the (unrealistic) assumption is made that the qubits 
are divided into physical blocks, each of length N c , which do not interact with each other; 
that is, the interaction Hamiltonian is applied separately to each such block. The blocks 
are intended to represent the "logical" qubits in a computer encoded for error correction. 
When the code is concatenated, the decoupling assumption is applied at the lowest (deepest) 
level: For instance, for a code concatenated once, a logical qubit is made of iV c x iV c physical 
qubits, and it is assumed that the N c lowest-level blocks do not interact with each other. 
While unrealistic, it is expected that this approximation should not modify the results very 
much, based on the numerical calculations made using a similar "decoupling of registers" 
approximation in ref. 0. 

In the following Section, the action of this Hamiltonian on a register encoded for error 
correction is analyzed, and then in Section III the analysis is extended to a concatenated 
code. The goal is to derive an expression for the "failure probability" per qubit, that is to 
say, the probability that the evolution under the Hamiltonian (0) may lead to a state which 
cannot be recovered by the error-correction protocol. Section II shows a few general features 
of this failure probability, such as its dependence on the initial state considered, and the way 
some error amplitudes add coherently and some do not, depending on the code considered. 
These results pave the way for the more complicated situation discussed in Section III. It 
is to be noted that the failure probability is calculated here for arbitrarily large interaction 
strength, for up to two-level-deep concatenated codes. 



II. RESULTS FOR SIMPLE DISTANCE 3 CODES 

Distance 3 codes can correct any of the three basic one-qubit errors represented by the 
three Pauli matrices a x , a y , and a z . They will typically fail whenever two errors happen 
simultaneously in the same logical block. 

The Hamiltonian (0) causes two phase errors in each physical block every time it is 
applied, but the actual damage it does to an encoded state turns out to depend on the code 
used, and on the state it acts upon. Two examples are considered in this section. 

The evolution operator resulting from (pj) is, for a single logical block, 

N c -l 

U{t) = Y[ (costft + ia zn a z>n+1 sin 5t) (3) 

n=l 

since the Pauli matrices for different qubits commute. Thus, a simple perturbative estimate 
of the failure probability would be as follows: if it is assumed that each term of the form 
G zn°z,n+\ leads to a distinct (orthogonal) uncorrectible error, there are N c — 1 such terms in 
(§), and the sum of the squares of their amplitudes is 

P P ert(t)= (N c - 1) (cos 2 5t) Nc ~ l sin 2 5t 

^(N c -l)(5t) 2 (4) 

for sufficiently small interaction strength and time St. 
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As it turns out, this expression is not, in general, correct, even for small St, because of 
two assumptions: that all terms of the form cr Z n a z,n+i lead invariably to fatal (uncorrectible) 
errors, and that the different errors are orthogonal. This is explored in detail in what follows. 



A. The 7-qubit code 

Consider the 7-qubit code introduced by Steane |§: 

|0 L ) = -L(|0000000) + 10001111) + 10110011) 

+|0111100) + |1010101) + 11011010) 
+|1100110) + 11101001)) (5a) 

|1 L ) = -=(11111111) + 11110000) + 11001100) 
v8 

+|1000011) + |0101010) + |0100101) 

+ |0011001) + 10010110)) (5b) 

Because the code is non-degenerate, all 7 states of the form a zn \0L) (with n = 1,...,7) 
are orthogonal to each other and to \0l)- Altogether, then, these eight states form an 
orthogonal basis of the space spanned by the eight kets appearing in (|5a|). But then, since 
any combination of products of a zn acting on |0^) does nothing but to change the relative 
signs in the superposition of these eight kets, it follows that it can always be written as a 
superposition of "one-error" states: 

7 

U(t)\0 L ) = C \0 L ) + Y l <rzn\0 L ) (6) 

n=l 

and this is always correctible. A similar reasoning applies to the state |1_l), which means 
that, if the initial state of the system is either \0l) or the evolution operator (^) does 
not introduce any fatal errors at all. 

On the other hand, the initial state of the system is much more likely to be a coherent 
superposition of \0l) and and for such states the situation is different, because, although 
error correction will restore \0l) and separately, it will also introduce a phase difference 
between them. The following useful relations can be easily established from the explicit 
forms (H): 

Oz\O z i\§l) = o^o+sIOlH -cr z3 |0i) 

O"z2O"23|0l) = (?zA(?zb\QL) = O zG 27 |0 l) '= ~ O"zl|0i) 

ct Z 3(Tz4\0l)= -v z7 \0 l ) (7a) 

CziO^lU) = ^z^ zQ \l L )= cr zZ \l L ) 
o ~z20 *3|1l) = o - z4 o - z5 \1l) = o - zG o z i\Il)= Czi|1l) 

crz3<7z4\lL)= cr z7 \l L ) (7b) 
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and it is easily seen from this that when error correction is applied to the state 
&zn <T z,n+i( a \0L) + the result will be the state (— a\0i) + /3| lx,)) , that is, a state 

with the wrong relative phase. Equations ([?]) also show that the assumption that different 
Czn^.n+i lead to orthogonal error states is incorrect: for instance, <7 2 i<7 22 and cr 25 <7 26 lead to 
the same error state. 

Using Eqs. (0) it is straightforward to derive the correct form of the approximate failure 
probability for small St. When only the terms of first order in St are kept in the expression 
(§), acting on a general state a|0i) + (3\1l), one obtains 

St (3a 2l + 2a z3 + a z7 )(-a\0 L ) + (3\1 L )) (8) 

which means that the probability that after error correction one might end up with the 
state -a\0 L ) + (3\1 L ) is (3 2 + 2 2 + l 2 )(5t) 2 = U(St) 2 . One might directly call this the 
failure probability (since the state — a|0^) + is generally not equal to the initial state 

«|0l)+/3|1l)) or define the failure probability as including explicitly the square of the overlap 
between the initial and the final state (see Eq. (|T0|) , below, for the precise definition), in 
which case one finds 

P ^ U(St) 2 (l - (\a\ 2 - \f3\ 2 )^ (9) 

where the ~ sign is a reminder that this is only valid for small St. The failure probability @ 
is largest for the symmetric linear superpositions (|0^) ± \1l))/\/2. Note how it is different 
from the result (§) (for N c = 7), obtained under the assumption that the errors generated 
by different <J zn a Ztn+ i were all orthogonal, and it is also not equal to the same expression 
with N c — 1 replaced by (iV c — l) 2 , which would correspond to a "worst case scenario" in 
which all the different errors added coherently. The coefficient 14 falls instead in between 
the two extremes of 6 (N c - 1) and 36 {(N c - l) 2 ). 

The definition of the failure probability P used above is 

p = i- |(^ol^(t)l*o)l 2 - E K^oK„t/(t)|^o>l 2 (io) 

n=l 

that is, one minus the total probability that, after a time t, the system may be found in 
either the original state l^o); or one of the N c correctible states o" 2 n |\l/ ) (there are, of course, 
other correctible states, of the form <7 x , n \^o) and a y:n \^/o), but they cannot be accessed by 
this interaction operator). 

Using this definition it is relatively straightforward to calculate P exactly for arbitrarily 
large St. The following relation is useful: 

(0 L Wz,n 1 ,---,Vz,n k U\0 L ) = (-l) k (l L \a z>ni ,...,a Z:nk U\l L ) (11a) 

(0 L \a z , ni ,...,a Zjnk U\l L ) =0 (lib) 

where k — 0, ... ,7. This follows directly from the fact that U involves always the product 
of an even number of o z operators, and that the code words used in |0^) are the comple- 
mentaries of the ones used in |1^). 
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For future use it is convenient to define 



/o = |(0 L |f/|0 L )| 2 = |(l i |[/|l i )| 2 (12a) 



7 7 

|2 l /i \tt\i \ |2 



a = e kolk^icui 2 = x: i(uic/iu)r (12b) 

n=l n=l 

Direct calculation then yields 

/„ = 1 - f x = cos 2 (5t) cos 2 (25t) cos 2 (35t) (13) 

and 

P(*)=/ 1 (l-(N a -| ) 9| 3 ) 2 ) (14) 

and it is easy to verify that (0) is the lowest-order term in the expansion of (fH| ) in powers 
of St. 

The failure probability (|T4] ) for the case \a\ = is plotted in Figure 1 (solid line) as a 
function of St, from <5t = to 7r/2. From to 7r, the curve is symmetric about 7r/2, as is to 
be expected from the form of the evolution operator (||). 

B. The 5-qubit code 

In 1996 Laflamme et al. |J introduced a 5-qubit code which can correct general one-qubit 
errors. A possible encoding is 

\0 L ) = -^=(-|00000) + |01111) - |10011) + |11100) 
v8 

+|10101) + |11010) + |00110) + |01001)) (15a) 

\1 L ) = -=(-111111) + |10000) + |01100) - |00011) 
v8 

-101010) - |00101) + |11001) + |10110)) (15b) 

The codewords can be obtained from the 7-qubit code by deleting any two qubits (in this 
case, the first two), but the pattern of signs is nontrivial. 

Note that the logical and 1 states are still superpositions of eight basis vectors each 
(where the basis vectors are simultaneous eigenstates of all the <J z ,n)', hence, the ten one-error 
states cr^nlOi) and cr Ztn \li), plus the two no-error states \0l) and |1_l) are not enough to 
span the space resulting from all the possible combinations of the 16 basis kets. As a result, 
there are now no states which exhibit fully correctible evolution (unlike the \0l) and |l^) 
did for the 7-qubit code). However, we still have the phenomenon of non-orthogonal errors. 
In particular, we find that 

<^2O" z3 |0l) = a z4 cr z5 \0 L ) (16a) 
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Cz2Cz3|1l) = 0" 2 40" 2 5|1l) 



(16b) 



and neither of these expressions is equivalent to a linear combination of single-qubit errors 
acting on \0 L ) or 1 1^) ; in fact, they are orthogonal to all such states, and to the no-error 
states \0l) and |1l) as well. 

The calculation of P{t), as given by the definition ([II]) is a little more involved, although 
the results ( |TTD still hold. Defining as before 

/o = |(OL|f/|0 L )| 2 = |(U|f/|U| 2 (17a) 

5 5 

/i = £ I(Ol|^|o l )| 2 = £ \(i L \u\i L )\ 2 (i7b) 

n=l n=l 

one finds, thanks to Eqs. (|n|), that P(t) can still be written as 

P(t) = l-f -f 1 (\a\ 2 -\[3\ 2 ) 2 (18) 
but this time, unlike for the 7-qubit code, one does not have / = 1 — A- Instead, 

fo = cos 4 (5t) cos 2 (25t) (19a) 



/i = sm 2 (6t) (1 + cos 2 {St) cos 2 (25t) 



(19b) 



The failure probability ([TJ) for the case \a\ = \(3\ is plotted in Figure 2 (solid line) as a 
function of 5t, from 5t = to tt/2. For small values of St Eq. (|T8"|) reduces to 



P ~ (<5i) s 



6-2 (la 



(20) 



The reason / ^ 1 - /i for the 5-qubit code is that, as mentioned above, the states |0l) 
and <J Zt n\0L,) do not form a complete basis of the space spanned by the eight kets which 
make up |0l). For the 7-qubit code, instead, they do, and, since the operator U acts entirely 
within that space (when the initial state is |0i)), any part of U\0l) which is not accounted 
for by the original state |0x) has to show up in one of the cx 2jn |0i,) (and the same for the 
initial state |1l)). 



III. CONCATENATED CODES 

A relatively straightforward way to obtain a code that can correct more than one-qubit 
errors is to concatenate a one-error correcting code, in principle, as many times (or levels 
deep) as necessary In a code concatenated once, each logical qubit consits of N c blocks 
of N c (lower-level) qubits each. This code can protect against general errors in any three 
lower-level qubits, and also against many other combinations of multiple lower-level errors. 
The minimum number of lower-level errors required to cause the code to fail is four, two in 
one upper level block (causing that whole block to fail) and two in a different one. 
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In this section I look first at the 7-qubit and 5-qubit codes concatenated once. A physical 
(lowest-level) qubit will be identified by two indices, n and m, running from 1 to N c : the 
first index denotes the upper- level block it belongs to, and the second index the qubit 's 
position within that block. With the assumption that different blocks do not interact, the 
time-evolution operator @ becomes 

N c N c -1 

U(t) = Yl II {cos 5t + ia nm a n , m+1 sin 5t) (21) 

n=l m=l 

From now on, for simplicity, the subscript z on a z will be dropped, as only a z is used. 



A. The 7-qubit code 

If all the errors generated by the individual terms in ([H]) were orthogonal, it would be 
a relatively simple matter to calculate the probability of an uncorrectible error to lowest 
order. It would be given by the terms of the form o" nimi cr„ 1)jni+ icr„ 2m2 (j„ 2)m2+ i, with ri\ ^ n 2 
and m 1; m 2 arbitrary. There are C(N C , 2) ways to choose the indices ni and n^, and for each 
choice there are (N c — l) 2 different ways to choose mi, rri2, so there is a total of 

(22) 




terms in fl2ll) which contribute to the failure probability in lowest order, and all these terms 
are proportional to (8i) 4 (after squaring the probability amplitudes). 

However, as has been shown already in the previous section, not all the terms enumerated 
in ( |22"D are mutually orthogonal. For the 7-qubit code, according to (0), there are two 
sets of values for mi (one set of three values, and one set of two values) which result in 
identical states. When the equivalent choices for m 2 are factored in and added, and the 
probability amplitudes of the corresponding states are squared, one finds, instead of the 
factor (N c - l) 2 = 36 in @, a factor 

(3 2 + 2 2 + 1 2 )(3 2 + 2 2 + l 2 ) = 14 2 = 196 (23) 

As a result, to lowest order, the failure probability for this two-level code is given by 

, 2^ 



P^196Q(^ 4 (l-(|c| 2 -|/3| 2 ) : 



4116 (<5t) 4 ( 1 - (H 2 - \f3\ 2 ) 2 ) (24) 



for a state «|0l) + (3\1l)- The dependence on a and f3 exhibited by (|24[) may not be 
not immediately obvious at this point (it will be derived later), but it shows that, for the 
concatenated code also, the individual states |0^) and |1^) are fully correctible. 

Note that the maximum error probability P^l x for the once-concatenated code, according 
to (p4[) , is related to the corresponding one for the plain code, (|9|), by the formula 



p(2) = ') (p(i) Y (25 ) 

max I r> I V max I V "/ 




s 



and the right-hand side of (|25|) is simply the probability that two blocks might fail inde- 
pendently, times the number of ways to choose two blocks. This is because, to this order 
(recall that these are only approximate results for small St), errors across different blocks 
are orthogonal, ie., they add incoherently. This can also be seen from Equation (|23|) : within 
each block there is some coherent addition of error amplitudes (which results in the terms 
3 2 and 2 2 ), but then the total error probabilities for each block are simply multiplied. 

This is an important result which generalizes to higher orders, as discussed later in this 
Section (subsection C). The key point is that all the errors which would cause a particular 
block to fail are orthogonal to all the errors that would cause a different block to fail, to 
lowest order. This is because to lowest order all it takes is two errors in one block to cause 
that block to fail, and because the original code is nondegenerate, which means that all 
expectation values of the form {^,i\o'iO'j\rii J ) (where n, £ = 0,1) within the same block are 
zero. 

Equation (|2"5j ) does not, however, hold for arbitrarily large 5t, when multiple (more than 
two) errors within the same block need to be taken into account. For a code concatenated 
only once it is still possible, albeit very cumbersome, to derive this failure probability for 
arbitrary 5t. The definition (pi]) has to be changed to 

P = 1 - |(* |t/(i)|*o>| 2 - Y.\^o\A n U(t)\V )\ 2 (26) 

n 

where the {A n } represent a maximal set of correctible error operators leading to orthogonal 
states. The difficulty is that now, unlike in the previous Section, not all possible correctible 
errors lead to orthogonal states. For instance, any number of errors within one single block 
are now correctible, but certainly not all of these lead to orthogonal states. Thus, even 
though we started out with a nondegenerate code, the concatenated code is degenerate [jlO 



As far as I can tell, the set {A n } in ( PB"[ ) must be constructed by inspection. For the 7- 
qubit code, concatenated once, I find that the maximum number of correctible errors acting 
on a single block and leading to orthogonal states is 15, which can be chosen to be 

O-nm (m = 1, • • • , 7) 

Oniow (m = 2, . . . , 7) 

0~n20~n3 

(27) 

Of these, the first 7 are, of course, fully correctible at the lower level, whereas the others, 
upon correction at the lower level, lead to what look like "block" phase errors at the higher 
level. 

The single-block errors (^7|) can be combined in a large number of ways among the 
different blocks to yield other, orthogonal, correctible errors. Fortunately, however, the 
actual number of such combinations which make a nonvanishing contribution to (^6|) is 
relatively small. To see this, divide the error operators acting on a single block into those 
involving an even number of a's (denoted as E n if they act on the n-th block) and those 
involving an odd number of ex's (denoted as O n ). Then, for a general multiblock error of the 
form A = E ni . . . E nk O mi . . . O mk , Eqs. ([II]) for a single block imply that the expectation 
value 
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(0 L \O ni ...O rik E mi ...E mi U\0 L ) (28) 

will only be nonzero if the corresponding combination (0L|cr ni . . . cr n jOz,) is nonzero for the 
single, non-concatenated, code. The same is true for the state and there are no cross- 
terms between |0l) and (also as a direct consequence of the single-block Eqs. (|Tl"])). 
Thus, to evaluate (|26"D , it is sufficient to identify the combinations of numbers {n\, . . . , n^} 
which lead to a nonzero (0^|(7 ni . . . c^jO^) for the 7-qubit code, and put odd errors, chosen 
from the list (p7|) , in the blocks n\, . . . , Uk, augmented perhaps with even errors, always 
chosen from the list (P7|), in any other blocks, provided, of course, that the total error be 
correctible. This quickly reduces the total of different possibilities to a manageable number: 
for instance, since an even error in one block would cause that block to fail, one cannot have 
even errors in more than one block at a time. 

For reference, for the 7-qubit code, the sets of values of 7ii,...,7i& to be used are: 
(1,2,3), (1,4,5), (1,6,7), (2,4,6), (2,5,7), (3,4,7), (3,5,6), the complementary set 
((4,5,6,7), (2,3,6,7), and so on), and the set (1,2,3,4,5,6,7). Then, possible error 
operators A n to be used in fl26|) might be products such as such as cr ljmi cr 2 ,m2 c,r 3,m3, or 
o"i,miO"2,m 2 o"3,m 3 o"4,i '4, m 4 ) or ^n,i^n,m- In terms of the functions f , fi defined in flTJ), the 
term(s) T, mi ,m2,m 3 I (Oiki.mi^ma^ma ? contribute a term f£ff to the expression (gBJ): 
the term T, mi ,m 2 ,m 3 ,m 4 |(0L|cri, mi cr2,m 2 c r 3,m3O"4,ic r 4,m 4 |0L)| 2 contributes a term f$f£ (since the 
product of two cr's in the same block is equivalent, except for a sign, to a single a acting on 
the appropriate state, cf. Eqs. (^)), and so on. 

Collecting all the terms, and keeping track carefully of which ones have opposite signs 
in the states |0^) and the states 1 1^) , one finds 

P {2 \t) = (l " fo - 7/oVi - 28/ 4 A 3 - 7 ft f? ~ 21/ 2 A 5 ) (l - (l«| 2 - \P\ 2 ) 2 ) (29) 



with /o and f% given by (0), which does reduce to ( p4|) for small St. 

The exact result ( p9|) is plotted as a function of St in Figure 1 (dashed line), for |a| = 
It is, perhaps, remarkable that it does not look all that different from the single-encoding 



result (14), except in the very small St region. It is also clear that Eq. (|25f) does not hold for 
large St, which means that for large St we do see the effects of errors across different blocks 
adding coherently. 



B. The 5-qubit code 

For the 5-qubit code, one again has to identify the appropriate maximal set of error 
operators leading to orthogonal errors. For a single block, one finds the following set of 15: 





(m = 


1,. 


.,5) 


&nl@nin 


(m = 


2,. 


.,5) 




(m = 


3,. 


-,5) 




(m = 


3,. 


-,5) 



(30) 



As in the previous subsection, it is found that only certain combinations of odd errors 
across different blocks lead to nonzero contributions to (|26"|) . These turn out to be much 
fewer than before, being limited to n 1 , . . . , n k = (1, 2, 5), (1, 3, 4), and (2, 3, 4, 5). 
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The result, valid for arbitrarily large St, turns out to be now 



P^(t) = l-(/ 5 + 5/ 4 (l - f ) + ft + 4/oA 3 (l - /i)) 

-(4/o/r + 6/ 2 A 2 -8/ 2 /r)(|c| 2 -|/?| 2 ) 2 (31) 



with /o and f\ given by Eqs. (|T9|). The factors of 1 — fo appearing in (|31| ) come from the sum 
of the single-block expectation values squared of (J\a m , m — 2, ... ,5 and cr 2 cr m , m — 3, ... ,5; 
these seven error operators, acting on the singly encoded \0l) or generate the seven 
orthogonal states necessary to complete a basis of the space spanned by the 8 kets appearing 
in (p~5|) . The factors of 1 — /i come from the sum of the single-block expectation values 
squared of o\o%o m , m — 3, ... ,5: these 3 error operators turn out to be the complement of 
the 5 ones appearing in the definition ( |16b| ) of f\. 
The limit of (|3lT) for small St is 



P (2 \t) ~ (St) 4 



360 - 24 la 



12 1/312 



(32) 



As usual, the failure probability is largest for the symmetric superposition state (|0l) ± 
|1l))/v / 2- The corresponding coefficient, 360, can be obtained also from the simple argu- 
ments used to derive the result (|24| ) in the previous subsection: it equals C(5, 2)(2 2 + l 2 + l 2 ) 2 , 
where the 2 2 stands for the amplitude square of the sum of the two error terms which lead 
to identical errors, as shown in Eq. ([TB]) (and the two l 2 stand for the other two orthogonal 
errors). 



The exact result ([Jl]) is plotted in Figure 2 (dashed line) as a function of St, for |a| = 
In this case, unlike for the 7-qubit code, the difference with the single-encoding formula (solid 
line, equation (|TBP) is substantial; the twice-encoded qubit has a substantially smaller failure 
probability for most values of St. The failure probability actually goes to zero at S = 7r/2, 
where Eq. fl2"T|) yields U = i Yin=i < J niO r n5- Upon inspection, it is found to be a peculiarity of 



the concatenated (two-levels deep) 5-qubit code that 

5 

1 1 On\Onh («|0l) + = oxxo^a-nozio^o^. (a\0 L ) + (3\1 L )) (33) 



n=l 



and the state on the right-hand side of ([J3|) is fully correctible, since it amounts to only one 
block failing. 



C. Results for arbitrarily deep concatenation 

Clearly, the exact calculation of P^ 1 ' for a code concatenated n levels deep, with n > 2, 
becomes much too cumbersome to be feasible. However, the results in the previous subsec- 
tions show that the leading term in St can be calculated using relatively simple arguments. 

Consider, for instance, a code concatenated three levels deep. It will fail if at least two 
of the high-level qubits fail, but these high-level qubits are themselves two-level deep codes. 
One can choose the two highest-level qubits that fail in any of C(N C , 2) ways, and for each 
of them the number of terms in the expansion of the evolution operator which lead to a 
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qubit failure are those computed for the previous level. So, for the 7-qubit code, we can 
immediately iterate (|24 ) to get 



,2 



4116(5t) 




4 >2 



(34) 

where the subscript "max" indicates that we are looking at the failure probability for the 
state that maximizes it (the linear combination (|0x) ±|1e,))/\/2), and Eq. (^) has been used 
for PW. 

Equation ([51) contains the essence of a recursion relation which yields the failure prob- 
ability at the n-th level of concatenation (n > 2) as 




(35) 

Figure 3 shows the approximate pW, calculated from (|35|), for n = 1, . . . , 4 (solid lines), and 



also the exact results (|T4| ) and (|29|), for comparison. Based on these examples, it would seem 
that Equation (35) overestimates the failure probability, which would make it a somewhat 
conservative estimate 

As encoding depth n increases, the failure probability (|35f) will, for a given 5t increase 
or decrease depending on whether C(7, tyP^lx is greater than or less than one. This yields 
the "threshold" value for (5t) 2 , below which virtually error-free operation can be achieved 
for sufficiently deep encoding. Using @, one finds the condition: 

{6t)2 < 21^14 = 3 ' 4 X 10_3 (36) 

(Figure 3 shows this as the point where all the curves cross.) This threshold (|36| ) is com- 
parable to some estimates for the correction of independent, random errors by concatenate 
codes, if one thinks of {St) 2 as a sort of failure probability for a single physical qubit. As 
pointed out in the Introduction, this is a natural interpretation (in spite of the fact that two 
physical qubits at a time are modified by this interaction), because of the result ([[]). 

A similar calculation yields the threshold result for the 5-qubit code (using fl20|) and 
C(5,2)=10): 

(5t) 2 < — *— = 1.7 x 1(T 2 (37) 
v ; 10-6 V ; 



IV. DISCUSSION 

The results in the previous section suggest that, in principle, the correction of even rather 
large undesired interactions between the qubits does not require special methods or resources 
beyond those needed to control independent, random errors arising from the interaction of 
the qubits with the environment. This observation is based, first, on the fact that estimates 
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of the threshold for fault-tolerant computation with independent environmental errors tend 
to be rather more restrictive than either fl36|) or (|37|) |||7|; and, second, on the fact that for 
most proposed quantum computing systems the unwanted interactions betwen qubits would 
easily satisfy the threshold conditions (|36|) or ([37|), provided that the time in between error 
corrections is not excessive. 

As an example, consider a system of nuclear spins, in a solid, for instance (such as in 
the proposal by Kane 0]), separated by ~ 150 A, and interacting via the magnetic dipole- 
dipole interaction, which would naturally lead to a Hamiltonian like (0). For such a system 
6 may be estimated as 

6 „ „ io-V 1 (38) 

Here /zat is the nuclear magneton. 

Without error correction 5t becomes already of the order of 1 after 100 s, and long 
calculations are impossible. But assume that one uses error correction and that complete 
error correction of an n-order encoded qubit requires a time t = c n T gate , where the time 
needed to perform an elementary logic gate for this system is about T gate = 10~ 5 s. Then 
the threshold is easily reached unless c n is exceptionally large, which would make error 
correction generally impossible anyway (a large c n would mean essentially that the computer 
architecture does not allow for a large degree of parallelism when correcting errors across 
different qubits). 

Of course, the results presented here have been obtained only for a specific kind of 
interaction Hamiltonian and a specific geometry, and would need to be modified for other 
cases (e.g., two-dimensional geometries, with more nearest-neighbors per qubit, or to include 
non- nearest neighbor interactions). There is also the simplifying assumption made from the 
outset, of neglecting interactions beween qubits in different blocks. These complications 
might require recalculations for specific physical systems and hardware configurations, but I 
do not expect them to change substantially the threshold estimates of the previous Section. 

For instance, for the kind of "linear" geometry considered here, including interactions 
between the last qubit of one block and the first one of the next block would not modify the 
failure probability to lowest order, since a single two-qubit error which straddles two blocks 
does not cause either block to fail; hence, to make any one block fail, more of these errors 
are required than of the kind of "internal" errors considered in this paper. 

The key result of Section III concerning the incoherent addition of failure probabilities 
for different blocks, to lowest order, should be quite general: it depends only on the base 
code being a nondegenerate distance 3 code and the interaction Hamiltonian being a sum 
of pairwise products of elementary error operators for individual qubits. 

The most important "real-life" complication which has been ignored here is, rather, the 
fact that, in practice, error correction will not be a clean and purely mathematical operation, 
but a physical one, involving interacting systems and therefore itself subject to error. The 
theory of fault-tolerant error correction has shown that even when this is accounted for it 
should still be possible to achieve "almost" perfect error correction, provided that the error 
rate per gate is below a certain threshold, but it is here where one finds the (possibly) 
very small error thresholds which constitute the main challenge for large-scale quantum 
computing. 
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FIGURES 



FIG. 1. The maximum failure probability for a qubit encoded with the plain 7-qubit code 
(solid line) and with the same code concatenated once (dashed line) as a function of the product 
of interaction strength and interaction time 

FIG. 2. The maximum failure probability for a qubit encoded with the 5-qubit code (solid line) 
and with the same code concatenated once (dashed line) as a function of the product of interaction 
strength and interaction time 

FIG. 3. Approximate results, valid for small St, for the maximum failure probability for a qubit 
encoded with the 7-qubit code concatenated 1 (no concatenation), 2, 3, and 4 levels deep (solid 
lines). The exact results (from Figure 1) for the cases 1 and 2 are shown as dashed lines. 
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